Trapped ions in laser fields: a benchmark for deformed quantum oscillators 
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Some properties of the non-linear coherent states (NCS), recognized by Vogel and de Matos Filho 
as dark states of a trapped ion, are extended to NCS on a circle, for which the Wigner functions 
are presented. These states are obtained by applying a suitable displacement operator Dh (a) to 
the vacuum state. The unity resolutions in terms of the projectors \a,h) (a, , |q,/i~^) (a, /i 
• are presented together with a measure allowing a resolution in terms of |q, h) (q, h\. Dh ict) is also 

I used for introducing the probability distribution funtion pA,h (2) while the existence of a measure is 

. exploited for extending the P-representation to these states. The weight of the n-th Fock state of 

the NCS relative to a trapped ion with Lamb-Dicke parameter 77, oscillates so wildly as n grows up 
?— j , to infinity that the normalized NCS fill the open circle 77"^ in the complex a-plane. In addition this 

prevents the existence of a measure including normalizable states only. This difficulty is overcome 
by introducing a family of deformations which are rational functions of n, each of them admitting 
a measure. By increasing the degree of these rational approximations the deformation of a trapped 
ion can be approximated with any degree of accuracy and the formalism of the P-representation can 
be applied. 
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The theory of certain one-parameter (q-) deformations of Lie algebras, the so called quantum groups has been of 
great interest in the last decade in several areas of physics. In 1989 Biedenharn [1] and McFarlane [2] independently 
, defined the q-analogue coherent state of a deformed q-oscillator, for which Nelson et al. [3] were able to obtain the 
' resolution of unity. Since then the properties of a class of deformations of the harmonic oscillator were considered by 
several authors (see f.i. [4]). Deformed quantum oscillators are represented by dynamical variables A, and Na 
satisfying the commutation relations, [yl,7V^] = A^ [A\Na\ = —A^ and [A, A'''] = f {Na), with f {Na) an arbitrary 
I ' real function of Na- All such variables are constructed in terms of single-mode field operators]^ a, and a^a. 
^ : In 1993 Crnugelj et al. [5] observed that the multiphoton interaction of a single mode laser field with a two level atom 
' is described by deformed-oscillator creation and annihilation operators which in combination with the pseudo-spin 
^ , atomic operators cr+ and a- , form the potential 



>< 



Wjc = A^a- + Aa+ 



used in the Jaynes-Cummings model (JCM) [6]. 



Jh In that period the J-C model was at center of the attention for the study of laser cooling of ions placed in parabolic 

. . ! traps, with the quantized center-of-mass motion of the ion playing the role of the boson mode, coupled via the laser 
to the internal degrees of freedom. In the case of cooling the operator A is represented by a combination of some 
power of the annihilation operator a times a function of n. When some bosons of the oscillator mode are destroyed 
the ion is excited to the upper level from where it decays radiatively. Cooling was investigated in Lamb-Dicke [7] and 
strong-sideband [8] limits, that is for ion excursions small compared with the radiation wavelength. This study led 
to the discovery of many intriguing effects connected with the nonclassical properties of the field, like as a long-time 
sensitivity to the statistical properties of the radiation field [9]. For example, the mean excitation number of the 
quantized oscillations of a ion driven by a squeezed field exhibited periodic collapses and revivals [10]. 

The interest for the vibrational motion of trapped ions was also motivated by the connection between the state of 
motion and the properties of the fluorescence spectra [11,12]. This link led some experimentalists to look for new 
non classical radiation states generated by trapped ions forced into some unusual vibrational states. In analogy to 
the preparation of nonclassical states of light in quantum optics several authors examined the preparation of the 



^wherever possible operators will be indicated by simple letters, except for the addition of a caret when confusion could arise 
with c-number quantities 
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ccnter-of-mass motion in a quantum state having no classical counterpart. Worthy examples were those of Cirac 
et al. [13] who considered the possibility of generating squeezed states of the vibrational motion by irradiating the 
trapped ion with two standing-wave light fields of different frequencies and locating the center of the trap potential 
at a common node of both waves. In all these cases the nonlinear dependence of A on a, and n, stemmed from the 
ion motion in the trap potential. 

de Matos Filho and Vogel [14] observed in 1993 that the center-of-mass state of a trapped ion driven by a two-mode 
laser field decays toward a dark state coincident with a nonlinear coherent state (hereinafter called NCS) of a deformed 
oscillator. This result brought new fuel to the study of deformed oscillators describing different classes of states arising 
in the trapped ion motion under the action of two or three fields detuned by multiples of the vibrational frequency 
(see f.i. [15] for nonlinear cat states). In the wake of this interest attention was paid to theoretical models of deformed 
oscillators, like those connected with excited coherent states and binomial states [16] . All these non-linear oscillators 
differ for the deformation function h{h) connecting the annihilation operator a to the deformed one A = ah{n). 

The ancestor of these realizations were the q-oscillators characterized by a deformation hq (n) = \f^^^^^ increasing 
with n. Contrarily the trapped ion deformation is a very irregular function of n, taking positive and negative values. 
What is worse, for some combinations of the Lamb-Dicke parameter t]"^ and n it can vanish or become infinite. As 
a consequence it is hard to capitalize on the work done for the q-oscillator for studying the NCS of a trapped ion. 
In particular, while for the q-case it has been found a measure resolving the unity, the same is not exactly true for 
the ion case. As a consequence the formalism of the Bargmann spaces [17], which has been extended from the linear 
oscillators to the q-ones, cannot be applied exactly to the ion case. In fact, it will be shown in the following that 
this can be done by considering a class of rational deformations which approximate to any degree of accuracy the ion 
deformation. In most experimental cases the statistical state of a trapped ion is limited to a finite number of Fock 
states so that these rational deformations may adequately approximate the ion deformation. Only in this "weak" 
sense it is possible to construct a ion-analogue of a Bargmann space, on which the deformed creation and annihilation 
operators are represented as multiplication by z and differentiation with respect to z, respectively. 

This paper is dedicated to an extension of the theory of the usual coherent states to NCS using as examples the 
deformation relative to the dark states of trapped ions. We start with a single-mode excitation field A (Sec. II), by 
discussing some properties of NCS, and introducing a deformed version (a) of the displacement operator (Sec. III). 
In Sec. IV we discuss some aspects of the resolution of unity for these NCS. The operator Dh (a) is used in Sec. V 
for associating the density matrix operator p to a linear functional pA,h (z) mapping the test function exp (az* — a* z) 
into the expectation value {Dh (cc)), by extending the construction of the antinormal probability distribution function 
[18]. The connection with the P-representation is also briefly examined. Section VI is dedicated to NCS on a circle, 
for which the Wigner functions are presented. Finally, the last section is dedicated to the dark states, arising when 
a trapped ion is driven by a bichromatic laser field. An asymptotic expression of the deformation and the relative 
factorial is obtained and its implication on the convergence of the NCS series is discussed. It comes out that it 
converges only for a in a circle of radius equal to the inverse of rj. On the other hand the weight of each Fock state 
can take values so large to prevent the resolution of unity in terms of normalized NCS. Some approximate expressions 
of the deformation are discussed together with the possibility of using these NCS for representing the ion statistical 
state. 



II. MOTION OF A TRAPPED AND LASER-DRIVEN ION 

We consider an ideal two-level ion of mass M constrained to move in a 3D harmonic potential. Taking the principal 
trap (x-axis) axis to coincide with the direction of propagation of the driving field, one quantum number suffices to 
label the vibrational states of the trap. The other two are traced out by summing over the corresponding degrees of 
freedom. 

The ion's internal and external degrees of freedom are coupled together by a light field ge.'^'^Lt+iv{t) periodically 
modulated at the frequency u of the ion trap 

E (x, t) = £e'"'^*+*'^(*\g {t) f (x) + h.c. 

where g {t) = 9 {t + ^) is a generally complex periodic function of frequency and h.c. stays for the Hermitian 
conjugate. The function / (x) stands for e"'*^^^ or sin (fc^x + </>) respectively for a progressive or standing wave, with 
the phase cf) determining the position of the trap potential with respect to the standing wave. 
We will dwell on monochromatic 

it) = e-^^^^'>' 
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and bichromatic driving fields 



9n+i — ^ ajv+1 

with the parameter N taking non-negative integer values, and aAr+i a complex coefficient depending on the amplitudes 
of the two waves. 

Now, introducing the Lamb-Dicke parameter rj = hkL j^lMfiv we put as usual e"'*^^ — e~"'(''"+'^'') . In the classical 
limit 77 is large and the absorption or emission of a photon will always cause some change in the vibrational state of 
the atom. In the non-classical Lamb-Dicke limit (LDL) of small 77, many photons may need to be absorbed or emitted 
before the atom changes vibrational state. For example in the sideband cooling experiment carried out by Diedrich 
et al. [19] the parameter r] was equal to 0.06. 

The Hamiltonian for a trapped ion interacting with a bichromatic field can be split in two parts 



where (h = 1) 



H — Ho + Hi. 



Ho = W12CT3 + '^fi (1) 



and, in the electric dipole approximation, 

H^nt = P {a-E* {x, t) + a+E {x, t)) 

When the Rabi frequency f?, relative to the laser induced transition between the ion ground and excited levels, is 
much smaller than the trapping potential frequency v, a perturbation expansion can be carried out in fi/z^, as discussed 
in Ref. [9]. This expansion allows a division into quickly and slowly varying density operator matrix elements, the 
former of which can be adiabatically eliminated. 

Arresting the calculation to the zeroth-order in fl/iy amounts to applying the rotating wave approximation. This 
approach can be easily pursued by switching to the interaction picture defined by the unitary operator Urw = 
exp [—i (i^Lfa -l- vfi) t] and retaining in the transformed hamiltonian H' the time-independent terms together with the 
slowly varying phase ip (t) of the laser field, 

H' ^ iA-ip{t))a3 + n{a_A + a+A^) (2) 

with A = UJ12 — llJl the detuning parameter, Q = e^^^^^pS the vibronic Rabi frequency and 

A = e-"'/2 g (t) / (e-'-'tat -t- e"'*a)] (3) 

the bar indicating the time average. 

Expanding the factor e^^^^'^ a'*+e a) power series in a and a^, introducing the operator 

- E i-^T - (4) 

with (n — m + 1)^^^ — (a^)™ a™ = n {n — 1) . . . {n ~ m + 1) and (77^) reducing in the Fock basis to the generalized 
Laguerre polynomials, we obtain respectively for progressive 



fc=0 

and standing waves 

sin [r; (at + a) + 0] = e~^"/^ ^ "^^T ('^ + ) b'' ^" + + ^" + («^)' e"'^'"*] (6) 

with k a positive integer and Cfc = ^ for fc = and ej, = 1 otherwise. 

For progressive (p) and stationary (s) monochromatic waves with g (t) — e^*'-^"'"^-'''* the operator A (see (||)) is 
given by 
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^N+l 

A« = Hr^in(<^+(A. + l)|)AW (7) 
while for two modes bichromatic driving fields 

^ (_,)A'+i si,, + (AT + /^^, (fi) a"+i + sin (<^) a^+ifo (n) . (8) 

III. NONLINEAR COHERENT STATES 

Coherent states were originally introduced as eigenstates of the annihilation operator for the harmonic oscillator 
[20]. They have been generalized (see [3, 4, 14, 16]) by labeling as nonlinear coherent states \a,h) the right-hand 
eigenstates 

A\a,h) = a\a,h) (9) 

of operator^ A of the form 

A = ah{h) (10) 

where h (n) is an operator-valued real function of the number operator. It is immediate to show that 

00 „ 

with [h {n)]\ = h (O) h (1) ■ ■ ■ h (n) and normalizing factor Nh^a 

1 



Eh ( |«|' 



expressed in terms of the entire function 



^o"!([M^)]!) 

referred to in the following as h-exponential in analogy with the q-exponential used in Ref . [3] . 

The deformation functions h (^n, ry^) associated to the dark states of the trapped ions are represented by the ratio of 
two Laguerre polynomials of argument equal to the Lamb-Dicke parameter 77^ so that they vanish or become infinite 
for some isolated combinations of 77^ and n. We are obliged to explicitly assume that this situation does not occur for 
the values of 77^ considered. 

In Sec. VII we will obtain an asymptotic expression of the weights of the Fock states occurring in the series 
expansion of the NCS relative to trapped ions. They take very large and very small values for increasing n, so that 
these NCS can be normalized only for a inside the circle 1/77. For convenience of discussion we shall ignore this 
problem by restricting our treatment here to normalized NCS states. 

It may be worth noting at this point that many of the foregoing formulas may be abbreviated by adopting a 
normalization different from the conventional one for the coherent state. If we introduce the symbol ||q!;/i) for the 
states normalized in the new way and define these as 



^for the sake of notational simplicity we will use the same symbol A for indicating fields of the form (^) and (^. 
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\\a-h) = N-X\a;h) 

and 

{a-h \\P;h) ^Eh (a*/3). 

Since the commutator \A\ A\ = nh^ [fi) — (n + l)h^{n + 1) is not a c-nmiiber it is worthy introducing the operator 
[16] 

4 - T^^' (12) 



With these alterations we have 



h{n) 



A \\a, h) = a \\a, h) , 
d 

n \\a, h) = a^T- \\a, h) 



da 
d_ 

da 



In addition 



||a, h) = (n) 4 ||a, h) = (^a^ ^ ||a, h) (13) 

In all the above r.h.s. the operators a, da and their combination are intended to act on the coefficients of the Fock 
states series. 

A. Displacement and deformation operators 

It is well known that coherent state \a) can be also introduced by displacing the Fock vacuum state |0) by means 
of the operator 

I?(a) = exp(-a*a + aa^) (14) 
due to its property of displacing the annihilation operator a by the generally complex quantity a, 

T) {a) oD {—a) ~ a — a 

Unfortunately, V {a) is unable to displace the deformed operator A. In alternative V [a) could be replaced by the 
unitary operator obtained by replacing in Eq. ( p^ a and by A and A^ respectively, but also this operator does 
not displace A by the complex quantity a. The difficulties in dealing with exponentials of linear combinations of A 
and A^ originate from the circumstance that their commutator is not a c-number. These problems can be overcome 
by using A^^ (see Eq. (p^)) in place of A^ and defining the "deformed" version of the displacement operator as 

Vh (a) = exp [-a* A + aA]}j = e^e-"*^e"^Ji = e-^e^^^^e""'-^ (15) 

{a) shares many properties of the standard operator D {a) as 

(a) = Vn (-a) 

and 



Vh {(i) Vh (a) = exp 
However, T^h (a) is not a unitary operator. 



^ {(ia* - P*a) 



Vh il3 + a) 



Vl (a) = VI' {a) = i-a) 
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so that it does not preserve the norm of a state. 

(a) and Vi (a) displace A and respectively by a and a*, 

h 

Vh (a) AVh (-a) = A-a 
Vi (a) A^Vi (-a) = A^ - a* (16) 

h h 

and Aj^ by a* 

Vh (a) AlVh i-a) - 4 - a* (17) 
Accordingly, the NCS \a; h) can be obtained by applying Vh {a) to the vacuum state, 

\\a-h)=e"< \Q)^e'-^Vh («) |0) 

In conclusion, the NCS ||a; h) can be obtained by deforming the usual coherent state ||a) by means of the deformation 
operator 

dh = Vh (a) V i-a) (18) 

namely, 

||a; h) = dh \\a) 

Although expressed as a product of operators depending on the complex parameter a, dh is independent of a. In a 
Fock basis it is diagonal with components equal to [h (n)]!^^. Since h (n) does not vanish, as already assumed, dh is 
not singular. 

Finally, we note that 

{m\Vh {a)\n)=^^^{m\V{a)\n) 

so that the matrix representation of V and Vh have the same diagonal part. 

A further remark is that the set of operators Vh (a) constitutes a Weyl system which does not lead to the canonical 
quantization for not being unitary. 



B. Nonlinear displaced Fock states 

In Sec. V we will use the Fock states displaced by Vh (a) (see Eq. (p7[)) 

\(pm,a,h) =Vh{a)\Tn) (19) 

which can be shown with the help of Eqs. ( |l6|) and (p^ ) to be the right eigenstates of the operator ^Aj^ — a*^ {A — a) = 
Vh (a) nVh {-a), 

(^Al - a*^ {A - a) |(/3m, a, h) ^ m \iprn, a, h) 

Analogously we can introduce the left eigenstates defined by 

{ipm, a, h\ (^Al - a*^ (A-a) = m (V'm, a, h\ 

which are obtained by displacing (m| by Vh (—a) , i.e. {ipm, a, h\ = (m| Vh {—a) . 

It is noteworthy that the left and right displaced Fock states are mutually orthogonal, 

{i;rn,a,h\ ipn,a,h) = 

for m ^ n. 

On the other hand these states can be also expressed in the form 
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T). ^ ' 



'm\ vm! 



(V'm, a, /i| = (m| (-a) = (a, /i| ..^f m?^^ = r, , ^ X! (?) " 

[/i(m)]!vm! [h {m)]W m\ \n J 



ml 



where \a, h, m) = Aj^ \a, h) and {a, h,m\ = {a, h\ stand for the deformed versions of the excited coherent states 
[20] (see also [21]). 

IV. RESOLUTION OF THE UNITY 

From the completeness relation of coherent states 

1 = ly |a) {a\d^a 



it descends 



Next, using the relation 



1 = — [ dh\a) {a\dh^(fa = — [ d^^ \a) {aldhd^a 
dj:' = V (a) V\ (a) = (Pi (a) V {-a)f = d\ 



the above resolution of unity can be expressed in terms of deformed coherent states 

^ = zl M~Zr W,h){a,h-'\d^a = ''7 \c^,h-'){(^Md''a (20) 



It goes without saying that this resolution holds true only if the NCS relative to the deformations h and 1 /h are both 
normalizable in the whole complex a-plane. 

For some deformations anyhow it is possible to obtain a resolution of unity in terms of projectors of deformed 
coherent states, i.e. to find a suitable element of measure d/j such that 



y ||a,/i)(a,/i||rf/i (21) 



d/i can be considered as an extension of the measure element dn = I"' d^a [17] for the linear oscillators. Since 
/ (to \\a,h) {a,h\\ n) dn must vanish for m 7^ n djjb can be put in the form 

where m/j {x) is a distribution satisfying the set of equations 

n\{[h{n)]\ f = j mh{x)x''dx (22) 

for every integer n. 

Treating n = s — 1 as a continuous variable the above relation represents a Mellin integral transform, 

gis)= fix)x'-'dx (23) 
Jo 

so that iTifi (x) is the Mellin antitransform of g (s) = T (s) {[h {s — l)]!)^- 

From the relation (/?, h\ \(i,h) = it descends that Ef^ (/3*a) is the self-reproducing kernel of the h-analogue 
of the Bargmann space [17], with respect to d^ 
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\Eh{f3*a)\^a"'dfi^ r 

In preparation of the discussion of Sec. VII it is worth remarking that replacing h by the deformation f3h the 
relative measure mph (x) is given by 



(24) 



This relation can be also used for expressing a thermal density matrix characterized by Boltzmann weight factors 
Pnn oc exp {—[3n) in the form 



ruh \a\ 



\a, h) (a, h\\ 



In Ref. [3] it was possible to obtain the resolution of unity for a q-oscillator by deforming both the derivative and 
the integral operators while a resolution for the so-called harmonious states was obtained in [23] . We will see in the 
following that for the trapped ion deformation the measure is a distributional Laplace antitransform which includes 
non-normalizable NCS. 

For a deformation approximated by a rational function of n, g (s) corresponds to the ratio of products of gamma 
functions. 



r (ai + s) • • • r + s) _ 
' (s) = — ^ — — r = r 



(a) + s 
{h) + s 



For A> B the relative antitransform is given by a combination of generalized hypergeometric functions 



(25) 



cPd 



(ci)„ ■ ■ ■ (cc)„ a:" ^ .nic+D+i) 



i-iy 



namely [24] 



ruh (x) 



A 
^=1 



(a)' - 
(6) - 



1 + - (a) 



(26) 



where (a)' — and (a)' — (a)' — — 1 stand for the sequences oi 
with the exclusion the ^-th term. 



, aA - an, and ai - a„ - 1, . . . , - flu - 1 



V. EXPANSION OF STATISTICAL STATES 



The same reasons that led [22] to express arbitrary states and operators in term of coherent states, suggest that we 
develop expansions in terms of NCS as well. 

Following [18] we introduce for a statistical state the deformed quantum linear functional 

Fh [a] = Tr {pVh {a)} = ^ Pmn {n\ ipm,Oi, h) (27) 

m 

In particular for a diagonal density matrix Fh [a] reduces to the standard F [a] . 
Using the unity resolution ( ^ ) Fh [a] may be rewritten as 

Fh[a]=e^ J exp (az* - a* z) ph.A (z) d^ z (28) 

where 

1 e^^^* 

Ph.A{z) = - .J \j Tr{p\z,h){z,h-^\} 

= -{z\ph\^) (29) 

TT 
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with 



Ph = ^pdh 



the deformed density operator. In other words Ph.A {z) stands for the generahzed distribution function of the deformed 
density matrix d^^pdh- 

For extending the definition of the characteristic functional F [a] — Tr {pD (a)} to a deformed oscillator we pay 
the penalty of loosing some properties of pA {z). In fact, ph,A {z) may take in general positive and negative values. It 
can be regarded as a generalized probability distribution function as long as the association between operators and 
functions is based on antinormal ordering, 

Tr [pGa (a, 4) } = y Ph.A {z) Ga {z, z*) cP z (30) 

In particular, for a diagonal density matrix pu.A {z) — Pa {z) while for p — \w, h) {w, h\ 

1 Eh (w*z) 
Ph,A [z) = -— — — exp [z {w - z)\ 
n bjh (w*w) 

Consequently, the transformation (E8|) applies if there exists the Fourier transform of exp [z* {w — z)] Eh {w* z). Anal- 
ogously, working with p = \w, h^^ | we arrive at the same conclusion for exp [z* {w — z)] Eh-i {w* z). This in 
turn implies that Eh-i (w* z) and Eh (w* z) cannot grow at infinity as quickly as exp (zz*). 



We recall that in the coherent states representation of a bounded operator O, the vanishing of 



O 



= in 



a domain of the complex plane of finite area implies the vanishing of O itself (see Refs. [18]). Since dh has been 

z), so that two deformed density matrices having 



assumed non singular the same theorem holds true for (^z Oh 

the same function ph.A (z) over some area of z, must coincide. In conclusion, Eq. (|29| ) establishes a one-to-one 
correspondence between the operator p and the function ph.A (z). 

When the deformation admits the unity resolution (^l|) a density matrix can be represented in several cases by a 
P-representation, 

Ph{a)\\a,h){a,h\\dp (31) 

in which Ph (a) can be regarded as a generalized probability distribution function as long as the association between 
operators and functions is based on normal ordering, 

Tr{pGN{A\A)} ^ f Ph{a)GN{a,a*)df, (32) 



When p is represented in the form ( pl| ) the master equation of p can be in many cases transformed in a master 
equation for Ph- This circumstance becomes particularly valuable in the study of the decay of an excited trapped ion 
toward the fundamental dark state. In this case we are faced for example with operators of the form 

A^P^ JP (a) mh (aa*) (^a-^^ \\a) {a\\ d'a 
I") ("II f-1 - ^4-) iP («) 



use having been made of Eq. (|13|) . Expanding P (a) mh (aa*) in power series of a™ (a*)" we see that 

oa \ oa J 
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VI. NONLINEAR COHERENT STATES ON A CIRCLE 



The above definition of NCS states (we will call them of order 1) can be extended to the eigenstates of the operators 
An+i of a more general form A^+i = a'^^^h (n) [25] and so the equation 

^jv+i \a, h,N + l,q) — a \a, h,N + 1, q) 

with TV > is considered. 

The eigenstate belonging to the eigenvalue a is A'' + 1-fold degenerate and q is an integer ranging from to iV. In 
terms of Fock states we have 

a\a,h,N + l,q) = Nc^haY^ = |^ (iV + 1) + g) (33) 

' •'t'oVinN+i)+qy.[hiiiN+i) + q)y: 

with the normalization factor 

1^ 1-2 ^ N 

UinN+l) + q)l{[h{l{N + l) + q)]\f 

where 

[h {l{N+l) + q)]l^ h{q)h{N + l + q)---h {I {N + 1) + q) . 

Such a state can also be expressed as a sum of NCS (see Eq.pl]). In facts, by introducing the function {n) 
defined recursively by 

a (N+l) + q)^ hiq-l)[h{l{N + l) + q)]l ^^^,^ 
^HiV+ij + gJ hiq)[h{l{N + l)+q-l)]r ^'^> 

we have also 

N 



h,N + l,q)=K^,^^Y. (34) 



k=0 



with e = exp (^j^^^ and t q ^ normalization coefficient. 

We have obtained that a NCS coherent state of order iV + 1 is decomposed in the sum of + 1 first order NCS 
of complex amplitudes a, ae'?, . . . ,a (e^)^ distributed uniformly on a circle. These states, referred to as "crystallized 
cats" in Ref. [26], were introduced for the linear oscillator [27] in the attempt to generalize the optical Schrodinger 
cats of harmonic oscillators. 

Using the deformed displacement operator we have also 



\a,h,N+\,q)=N' 



t N \ 

5] (£*)"'= (ae^) |0) 



In conclusion, the Hilbert space is the direct sum of + 1 spaces H = Tio (B Tii (B ■ • ■ (B Ti-N, (<Z = 0, 1, . . . iV), each 
of them having for basis the Fock states \£ {N + 1) + q), as / e (0, ...oo). 

For = 1 the fundamental states of Ho and T-Ci are respectively the even and odd Schrodinger cats. It will be 
shown in a following paper that when the radiative damping is negligible, the initial density matrix separate in the 
product of two matrices evolving respectively toward the even and odd Schrodinger cats. 

A. Wigner function 

The Wigner function [27] relative to these states on a circle can be shown to be given for a generic integer A^ by 
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w 

[/lA,+l (l iN+l) + q)]\ [hN+i {V {N + l)+ q)]\ [V {N + I) + q]\ 

2 



(35) 



q+ip 



Analogously for the Husimi-Kano [27] Q-function QAf+i,g(q,p) — 



with 



QN+i.q{q,p) = N^,h„+,,qe~ 



^^^J a coherent-state vector. 



a, 9, hN+i 

liN+l)+q 



E 

1=0 



[I iN+l) + q]\[hN+i {I (iV + l) + g)]! 



(36) 



For N = 1 these states reduce to even (q — 0) and odd {q = 1) Schrodinger cats [15]. In Ref. [28] is examined the 
squeezing and antibunching effects by using the function hi(n) introduced in [14] for the NCS. We will see in the 
following (see Eq. (^)) that the nonlinear cats representing the dark state of a trapped ion are properly described 
by the deformation h2 (n) — L^_2 (v^) /[n{n — 1) i„_2 (??^)]- 

In Fig. 1 we show the Wigner functions for nonlinear even Schrodinger cats of amplitude a = 3.5 (real) and different 
parameters r/. In the linear case (?7 = Fig. 1-a) the quantum interference is localized around the origin. The two 
coherent gaussian peaks are circularly shaped. 

For increasing rj the nonlinearity flattens the interference pattern while the central interference fringes, particularly 
their negative part, become more evident. This is essentially due to a reshaping of the coherent contribution peak 
from a gaussian-like nearly circular shape to an elliptical one with minor axis parallel to the direction connecting the 
two coherent peaks. These come closer to the origin and the region where the Wigner function is non-zero shrinks 
notably. 

A further increase in rj causes the progressive coming closer and closer of the main peaks, while the interference 
fringes become more localized. For higher rj there are some interference fringes spreading over the two coherent peaks. 
This phenomenon is dominant for very high 77 values (Fig. 1 b,c,d rj = 0.5) where the main peaks come into the 
interference region and the coherent character of the two states forming the cat is no more distinguishable. The 
interference area become larger than in the linear case and a circular symmetry of the interference pattern become 
evident. 

In Fig. 2 we present two NCS on a circle formed by the superposition of 3 and 4 NCS (a = 3.5 and 77 — 0.33). 



VII. DARK STATES 



Cirac et al [13] first proposed in 1993 a scheme for preparing coherent squeezed states of motion in an ion trap 
based on the multichromatic excitation of a trapped ion. Using two waves with beat frequency equal to twice the trap 
frequency a "dark resonance" appears in the fluorescence emitted by the ion, the ion is placed in a squeezed state. 
Similar "dark states" produced by a bichromatic field with beat frequency equal to the trap frequency were studied 
by Vogel et al. in 1996 [14] and identified as nonlinear coherent states. 

We will consider in the following a beat frequency which is a generic multiple of the trap frequency and the ion 
dark state is described by a generalized coherent state on a circle. We will consider a bichromatic field of the type 

Ab,N+i = Jn+i (n) a^^^ - aw+i/o (n) ■ (37) 
for which the dark state satisfies the equation 

that is (see Eq. (H)) 

\i'dark) = \aN+l,hN + l,N + l,q) 

with 
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hN+1 (n;r,2) = {N + l)l- 



tN+1 
^n-N- 



N) 



N+l 



(38) 



and 



no {N + iy. 

N+l {-i7])^+^ 



In short the dark state is the superposition of A'^ + 1 nonhnear coherent states which are equidistantly separated 
from each other along a circle with modulation factor e*"' = exp 
In particular 



/ 2-Kik 



0- 



hi (n;Tf) 
/i2 [fi; rl^) 



2Ll 



n{n-l) Ln-2 



(39) 



and 



.no 



a2 = -7^^ 



rio 2 



Expressing the Laguerre polynomials by their asymptotic expression 
depending on the product nr]^ only, 

TT 

4 



'nrfhi (n; 77^) tends, for n ^ 00, to a function 



'nrfhi {n; rj^^ = tan 



(40) 



Note the oscillating behavior of the eigenvalues E {n) ~ tan^ (^^/nrf — ^ /yy^. This circumstance implies that 

each eigenstate is encompassed by an infinite countable set of eigenstates of slightly different energies. For exceptional 
values of rj some eigenvalues can vanish. In these cases the series representation of the relative NCS looses its meaning. 
The behavior of E (n) as n ^ 00 has strong implication on the resolution of unity, as we will see in the following. 
As a consequence of (|40|), the logarithm of the factorial \h\ (n, ry^)] !n! times tends asymptotically to 



\og{[hl (n;ry2)]!n!r;2" 

n — ^00 1 



8r/2 



log 



tan 



dz + const 



const 



where u {x) is an oscillating entire function 

2 



u{x) =- {-Ax^ [Li2 {e~'^) + Li2 (-e*'^)] + [Li^ (e"*'^) - Li^ {-e'^)]) 



EM) 



k=0 



k Ax {2k + 1) - 1 
~(2k + Tf 



cos [{2k + 1) Ax] 



(41) 



with Lin {z) 
function of \ 



= YlkLi^^l^^ polylogarithm function. According to ( p| ) log ([/if (n, ry^)] In!??^") is an oscillating 
nrf , with the envelope expanding proportionally to a/ nrf , as shown in Fig. 3. This behavior is confirmed 
by the exact expressions of log ( \h\ (n, 77^)] Inliy^") plotted in Fig. 4 versus n for = 0.01, 0.02, 0.1 and 0.2. 

In particular there exists a countable infinite sequence of values of mf' for which \h\ (n, ry^)] \n\rf'^ is close to 

exp so that the NCS can be normalized only for |ary| less than one. In other words, contrarily to the linear 

coherent states the NCS relative to a trapped ion fill the open circle l/yy in the complex plane. As 77 — > the domain 
of existence tends to the whole complex plane, as for the linear coherent states. While these states are normalized 
for arj inside the unit circle, the scalar product of \ai,h) and |a2,/i) is defined even if one the numbers aior a2 
has an arbitrary modulus as long as the product of the moduli is less than l/r/^. A similar situation occurs for the 
harmonious states [22] described by the deformation function h {n) ~ \/ ^/n. 

For generic combinations of vibrational excitation levels and parameter if' the ion rovibronic dynamics fully displays 
its nonlinear character. An example of this feature has been seen above in connection with the discussion of the Wigner 
functions of some nonlinear Schrodinger cats and states on a circle of order 3 and 4. 
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A. Unity resolution 



Being these NCS restricted to values of a such that \arj\ < 1 the Melhn transform ( p3| ) reduces to the single-sided 
Laplace transform 



lim 



Then, g (s) is the right-sided Laplace transform of the distribution rrih (j] *) and tends asymptotically to the 



analytic function cxp 77 '^v (^\/sr^^ 



of s in the half-plane Re s > 1 and is bounded according to 
gis) <K e-^^"' 



with T a negative infinitesimally small constant and K a constant. Consequently rrih is a distribution with support 
bounded on the left at t = T < (see Ref. [30], corollary 8.4-la.). This means that it is necessary to include in the 
unity resolution un-normalizable states of amplitude \a\ > rj~^. 



B. Approximate deformations 



The difficulties in dealing with this deformation can be overcome by using approximate deformations. This is 
justified by the circumstance that in laser cooling experiments one deals with ions occupying a finite number of 
vibrational levels. 

In particular when the parameter rj^ is not very large the deformations {n) can be approximated by a few terms 
of the series expansion 



hi {n-ri^) = 1 + 



-V 



6 



48 



h2{n-if) =l + -{n-2)Tf + 



lln^-39n + 34 , 19n-^ - 96Ti^ + 159?i - 86 



24 



-■q 



60 



-7? 



(42) 



Approximating the deformation hi (n) by 1 + ^^^rj^ we have that ^ (x) coincides with the Mellin antitransform 



of F (s) + s — 1^ Before using Eq. (26) with A = 3, B = and ai = 0, a2 = = ^ — 1, we have to remove 



the degeneracy 02 — 03 by evaluating Eq.(g6|) for 02 — 03 = e and letting e 0. Since hm^^o T (e) + L (— e) 
— 2sinh(7£) /e = —2j with 7 = 0.57721 the Eulcr's constant, then 



ajTO/ij (a;) = [02] ; 1 - a2,l - 02; - — 

L ^2. 
At the same time the h-exponential reads 

Eh {v) = 0F3 



27r [-02] F [202] 0F2 



; 1, a2 + 1, 0,2 + 1;^ 



; 1, 1 -I- a; ^ 



Expanding the generalized factorial [hi (n, 77^)]! to the second order in 77^ the NCS \a,hi) can be expressed as a 
combination of excited coherent states 



|a, /ii) = (^1 + r,^ - j |a) + (^^ + ^ j a |a, I) + ^-J- + J-^ a' \a, 2) - -^a' \a, 3) + O {r,^) 

A better approximation can be obtained by representing the Laguerre polynomials L„_i (77^) by a finite sum of 
powers of rf, 



K 



fc=0 



J2k 



k k\ 
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so that hi {n-^rf^ (see Eq. (|39|)) can be replaced by a rational function 



I , Pa+1 - 1) - Pa+1 {n) „. n.^i [n + ai) 



7^^PB{n-l) 



(43) 



For using the characteristic function ph^A we should choose A = B, while for introducing the P-representation A> B. 
In particular, Cirac et al [8] [13] and Blockley et al [9] have expanded the exponential of Eq. (||) up to the second 
order in 77, their case corresponds to A = B = A. 

If the roots — and —hi are not integer we have (cf. Eq. (p5|)) 



(«) 



{a) + s 
{b) + s 



For this class of deformations the measure is given by the combination of the generalized hypergeometric functions of 
Eq. (p^), subject to the precaution of removing the degeneracy of the coefficients ai,bi. 

In particular for A = B = 1 the Husimi-Kano Q- function (36) relative to the state \z,h) reduces to 



Qi{w) a e-l"'-"l'/2 luj^J^^ {uj)f 



|ot|— >-oo 



1 



cot {VTV) 



with v = + 



1 = ^zw* and C, — (iff)^''- Consequently, the projector \z,h) {z,h\ is represented in terms of the 
undeformed coherent states by a P-representation containing derivatives of the Dirac function of the very high order 
V. This confirms the advantage of using the representation (p3l|). 

For a finite rank density matrix {x) can be represented by a finite combination of Laguerre polynomials 



Imposing the condition 



m/i (x) = e ^ ^ m„L„ {x) 
for Q <n < nmax yields 



(44) 



n=0 



{[h{m)W 



In Fig. 5 we have plotted these approximate measure functions for different values of {— 0.015, 0.0156, 0.0158, 0.016), 
representing density operators relative to ions excited up to the level n — 50. For these values of rf' the inclusion of 
a larger number of terms (n > 50) would lead to measures taking negative values. 



VIII. CONCLUSIONS 



The vibrational steady states of ions placed in a parabolic trap and driven by bichromatic fields detuned by multiples 
of the vibrational frequency provide a class of realizations of the nonlinear version of the so-called coherent states on a 
circle. The most well known example is that of the nonlinear Schrodinger cat states. As for the linear case also these 
states can be decomposed into finite sums of nonlinear coherent states, which can be considered as the building blocks 
of the vibrational wavefunctions of systems driven by laser fields detuned by multiples of the vibrational frequency. 

This class of states is well described by the Wigner function, which has been computed for states on a circle of 
degree two (cats), three and four. The relative patterns show a dramatic dependence on the Lamb-Dicke parameter 
rj, which measures the degree of departure from the linear case. This behavior is due to the irregular dependence of 
the ion deformation function on the Fock state index n. 

With the aim of investigating the possibility of extending some mathematical tools of the linear coherent state 
theory to NCS it has been introduced a deformed displacement operator (a), which in analogy with the linear 
one generates the NCS |a, /i) by displacing the Fock vacuum state (|a, /i) = Vh (a) |0)). The penalty paid for this 
extension is the loss of the unitarity. However, it allows the construction of a linear functional, which can be used for 
representing density operators by means of a generalized probability distribution function pA,h {z)- 

The peculiarities of NCS connected with trapped ions become evident when the deformation factorial is analyzed 
for very large n. The weight of the n-th Fock state contributing to a NCS exhibits an almost periodic behavior by 

taking very large and very small values of the order of e*'^"^/'"'^ . Consequently, the NCS can be normalized only for 
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a filling the open circle l/rj in the complex a-plane. In addition, this behavior prevents the existence of a regular 
measure for resolving the unity. These pathologies mark the difference with q-oscillators whose deformation function 
is an increasing function of n. 

For extending to these NCS the P-representation formalism it is necessary to replace the deformation with an 
approximate one, for which there exists a measure. These approximate NCS can be constructed in different ways. 
Two examples are provided. In the first case the deformation is represented by a rational function of the occupation 
number, obtained by truncating the Laguerre polynomials to some order in ry^. By increasing the degree of those 
rational functions it is possible to represent accurately the actual deformation for occupation numbers extending up 

to infinity. The respective measures are the Mellin antitransforms of Gamma function products T ^^^j and are given 

by combinations of generalized hypergeometric functions. In the other case the measure is represented by a finite 
combination of Laguerre polynomials. In this way it is possible to represent exactly the factorials up to a given level 
n, although it is not possible to obtain in general a positive definite measure. 

The NCS can provide a basis for studying the trapped ion evolution by representing the statistical expectation 

values of either antinormal Ga oi' normal Gn ^A, Aj^^ products of A and as integrals of the probability 

distributions ph,A (z) or (z) times the classical functions Ga (z, z*) and Gn (z, z*). Another possibility consists in 
transforming the density matrix master equation in an equivalent equation for the P-representation. This problem 
will be addressed in a more systematic way in a forthcoming paper. 

Before concluding, it is remarkable that the deformation used for the q-oscillators, the ancestors of the NCS, is based 
on the same transformation used by Heine [31] a century ago for generalizing the Gauss hypergeometric function. 



ACKNOWLEDGMENTS 



V.I.M. thanks the University of Naples "Federico 11" for the kind hospitality and the Russian Foundation for Basic 
Research for the partial support under Project No. 99-02-17753. 

References 

[I] L. Biedenharn, J. Phys. A22, L873 (1989) 

[2] A. MacFarlane, J. Phys. A22, 4581 (1989); see also C. P. Sun and H.-C Fu, J. Phys. A22, L983 (1989); M. 
Chaichian and P. Kulish, Phys. Lett. B232, 72 (1990) 

[3] R. W. Gray and C. A. Nelson, J. Phys. A23, L945 (1990); A. J. Bracken, D. S. McAnally, R. B. Zhang and M. 
D. Gould, J. Phys. A24, 1379 (1991); B. Jurco, Lett. Math. Phys. 21, 51 (1991); C. A. Nelson, "Novel imphcations 
of the q-analogue coherent states", in "Symmetries in Science VI", ed. B. Gruber, Plenum Press, N.Y., p. 563, (1993). 

[4] V. I. Man'ko, G. Marmo, E.C.G. Sudarshan, and F. Zaccaria, Physica Scripta 55, 528 (1997); V. I. Man'ko 
G.Marmo, and F. Zaccaria in "Symmetries in Science", Ed. B. Gruber and M. Ramck, Plenum Press, N. Y. (1997); 

[5] J. Crnugelj, M. Martinis, and V. Mikuta-Martinis, Phys. Lett. A 188, 347 (1994); Phys. Lett. B 318, 227 
(1993); Phys. Rev. A 50, 1785 (1994); 

[6] E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 (1963); 

[7] B. Buck and C. V. Sukumar, Phys. Lett. 81 A, 132 (1981); J. Phys. A 17, 885 (1984); V. Buzck, Phys. Rev. 
A 39, 3196 (1989), G. S. Agarwal, J. Opt. Soc. Am. B2, 480 (1985); C. C. Gerry, Phys. Rev. A 37, 2683 (1988); 
C. V. Sukumar and B. Buck, Phys. Lett. 83 A, 211 (1981); A. S. Shumovsky, Fam Le Kien and E. I. Aliskenderov, 
Phys. Lett. A 124, 351 (1987); 

[8] J. J. Cirac, R. Blatt, P. ZoUer, and W. D. Phillips, Phys. Rev. A 46, 2668 (1992) 

[9] C. A. Blockley and D. F. Walls, Phys. Rev. A 47, 2115 (1993); C. A. Blockey, D. F. Walls, and H. Risken, 
Europhys. Lett. 17, 509 (1992) 

[10] J. H. Eberly, N. B. Narozhny, and J. J. Sanchez-Mondragon, Phys. Rev. Lett. 44, 1323 (1980); N. B. Narozhny, 
J. J. Sanchez, and J. H. Eberly, Phys. Rev. A 23, 236 (1981); P. Meystre and M. S. Zubairy, Phys. Lett. 89 A, 390 
(1982); C. C. Gerry, Phys. Rev. A 37, 2683 (1988); J. R. Kukeinski and J. L. Madajczyk, Phys. Rev. A 37, 317 
(1988); 

[II] F. Dicdrich and H. Waher, Phys. Rev. Lett. 58, 203 (1987); M. Schubert, I. Siemers, R. Blatt, W. Neuhauser, 
and P. E. Toscheck, Phys. Rev. Lett. 68, 3016 (1992) 

[12] W. Vogel, J. Phys. B: At. Mol. Phys. 16, 4481 (1983); W. Vogel and Th. Ullmann, J. Opt. Soc. Am. B3, 

441 (1986) 

[13] J. J. Cirac, R. Blatt, A. S. Parkins, and P. Zoller, Phys. Rev. Lett. 70, 556 (1993); J. I. Cirac, P. Zoller, Phys. 
Rev. Lett. 74, 4091 (1995). 



15 



[14] R. L. de Matos Filho and W. Vogel, Phys. Rev. A49, 2812 (1994); Phys. Rev. A 54, 4560 (1996); 
[15] S. Mancini, Phys. Lett. A 233, 291 (1997) . 



[16] X-G Wang and H-C Fu, |quant-ph/9903013| vs. Nov (1999). 
[17] V. Bargmann, Commun. Pure and Appl. Math. 14, 187 (1961). 

[18] C. L. Mehta and E.C.G. Sudarshan, Phys. Rev. 138B, 274 (1965); see also J. R. Klauder, J. Math. Phys. 5, 
177 (1964) and T. F. Jordan, Phys. Lett. 11, 289 (1964) 

[19] F. Diedrich, J. C. Bergquist, Wayne M. Itano, and D. J. Wineland, Phys. Rev. Lett. 62, 403 (1989). 

[20] G.S. Agarwal and K. Tara, Phys. Rev. A 43, 492 (1991); V. V. Dodonov, Ya. A. Korennoy, and V. I. Man'ko 
and Y.A. Moukhin, Quantum and Semi classical Optics 8, 413 (1996); 



[21] S. Sivakumar, |quant-ph/9806061 



[22] R. J. Glauber, Phys. Rev. Lett., 10, 84 (1965); E. C. G. Sudarshan, Phys. Rev. Lett. 10, 277 (1965); 
[23] E. C. G. Sudarshan, Int. J. Theor. Phys. 32, 1069 (1993) 

[24 ] L. J. Slater, "Generalized Hypergeometric Functions", Cambridge University Press, Cambridge (1966); Proc. 
Camb. Phil. Soc. 51, 577 (1955) 

[25] Jinzuo Sun, Jisuo Wang, and Chuankui Wang, Phys. Rev. A 44, 3369 (1991); 
[26] O. Castanos, R. Lopez-Peiia, V.I. Man'ko, J.Russ. Laser Res., 16, 477 (1995) 

[27] J. Janszky, P. Domokos, and P. Adam, Phys. Rev. A 48, 2213 (1994); Ren-Shan Gong, Phys. Lett. A 233, 
297 (1997). 

[28] E. Wigner, Phys. Rev., 40, 749 (1932); K. Husimi, Proc. Phys. Math. Soc (1965). 

[29] B. Roy, Phys. Lett. A 249, 25 (1998); B.Roy and P. Roy, Phys. Lett. A 263, 48 (1999) 

[30] A. H. Zemanian, "Distribution Theory and Transform Analysis", McGraw-Hill, New York (1965). 

[31] E. Heine, Handuch die Kugelfunctionen, (1898) quoted by Slater (Ref.[24]). 

Figure captions 

Fig. 1: Top and side views of the Wigner function (Eq. |3^) relative to linear (a) and nonlinear (b,c,d) even 
Schrodinger cat states for different values of a = 3.5 (real) and two different values of rj. (a) 77 = (linear case), the 
two coherent peaks are almost circularly shaped the interference pattern is axially symmetric along the axis defined 
by the center of the coherent peaks; (b - top-view, c - side-view and d - air-view) rj = 0.5, the gaussian peaks are no 
more distinguishable from the circularly symmetric interference pattern. 



Fig. 2: Top view of the Wigner function (Eq. 35) for 3 and 4 NCS states sitting on the circle, a ~ 3.5 (real) and 
r] = 0.33. This value of the nonlinearity rj leads to a smooth reshaping of the coherent peaks from a nearly circular 
shape to an elliptic one. 

Fig. 3 Asymptotic expression of u{^/z) versus z (see Eq. (^)) 

Fig. 4 \og{[hl (n, r/2)] In!??^") versus n for 77^ = 0.01 (a),0.02 (b),0.1 (c) and 0.2 (d). 

Fig. 5 Measure rrih (x) versus x obtained as a combination of 50 Laguerre polynomials (Eq. (^)) and four different 
values of rj^ = 0.015, 0.0156, 0.0158, 0.016 
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